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Abstract 

The approach of nonrelativistic QCD (NRQCD) factorization was proposed to study inclusive pro- 
duction of a quarkonium. It is widely used and successful. However, a recent study of gluon fragmen- 
tation into a quarkonium at two-loop level shows that the factorization is broken. It is suggested that 
the color-octet NRQCD matrix elements should be modified by adding a gauge link to restore the fac- 
torization. The modified matrix elements may have extra soft-divergences at one-loop level which the 
unmodified can not have, and this can lead to a violation of the universality of these matrix elements. 
In this letter, we examine in detail the NRQCD factorization for inclusive quarkonium production in 
e + e~ annihilation at one-loop level. Our results show that the factorization can be made without the 
modification of NRQCD matrix elements and it can also be made for relativistic corrections. It turns 
out that the suggested gauge link will not lead to nonzero contributions to color-octet NRQCD matrix 
elements at one-loop level and at any order of v. Therefore the universality holds at least at one-loop 
level. 
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A quarkonium system provide a unique place to study the dynamics of QCD, because a quarkonium 
mainly consists of a heavy quark pair QQ and they move with a small velocity v. An extensive review 
of quarkonium physics can be found in pQ. A decade ago the approach of NRQCD factorization was 
proposed to study inclusive production of a quarkonium 2J . In this approach the production of a QQ pair 
can be studied with perturbative QCD because the mass m of Q provides a large scale, the formation of 
the QQ pair into a quarkonium is characterized with NRQCD matrix elements by an expansion in v. In 
order to have prediction power these matrix elements should be universal, i.e., they do not depend on 
how the QQ pair is produced. This leads to the NRQCD factorization. This approach is widely used 
and successful. Especially, it can systematically take higher Fock state components of a quarkonium, 
including those components which contain a heavy quark pair in color-octet, into account. A striking 
success of the approach is to explain the ip' anomaly at Tevatron [2] by taking color-octet components 
into account 

Although the approach is successful, a complete proof of the factorization does not exist. Studies of 
various processes at one- loop level really show that the factorization hold at one-loop level |51Rfl 171151 1§1 ITU]. 
But a recent study shows that the factorization for gluon fragmentation into a quarkonium at two-loop 
level is incomplete, indicated by that some uncancelled infrared(I.R.) divergences are not matched by 
NRQCD matrix elements. To restore the factorization, a gauge link is introduced to modify color-octet 
NRQCD matrix elements for matching these I.R. divergences|llj. At one-loop level, this gauge link can 
generate some extra I.R. divergences in the modified matrix elements than those unmodified. This can 
affect the factorization at one-loop level in the cases studied before and can lead to a violation of the 
universality of color-octet NRQCD matrix elements. 

It is the purpose of the letter to examine if the universality is lost and if the factorization can be 
done in inclusive production of a quarkonium through e + e~-annihilation through a color-octet QQ pair. 
Our analysis includes not only the leading contributions in the v-expansion, but also the contribution 
from relativistic corrections at order of v 2 . We show at one- loop level in detail how soft divergences are 
cancelled or matched by color-octet NRQCD matrix elements without the gauge link suggested in 
Our study also shows that the relativist correction for a color-octet QQ pair can also be factorized in 
the same way. To our knowledge, there is no known example studied at one-loop level to show that the 
NRQCD factorization holds. 

We consider the inclusive production of a quarkonium H in the process 

e + e - _, 7 * h + X, (1) 

where the virtual photon is with the momentum q 7 and q 2 is much larger than the square of the quarko- 
nium mass. For this process we need to calculate the tensor 

T^(P,q 7 ,H) =J2 f d 4 xe iq T x (0\J v (0)\H(P) + X)(H(P) +X\J"(x)\0) (2) 

where the quarkonium carries P and is the electric current. The NRQCD factorization in [2j suggests 
the tensor can be written in a factorized form: 

T^(P,q 7 ,H) = F^( 1 S )-(0\Os( 1 So,H)\0)+G^( 1 S )-(0\P s ( 1 S ,H)\0) 

+F^( 3 P ) • (0|O 8 ( 3 P ,#)|0) +F^ kl ( 3 P 2 ) ■ (0\oi ijh{kl} ( 3 P 2 ,H)\0) + • • • , (3) 

where the matrix elements in the right hand side are defined with NRQCD fields and can be found in 
We only consider the production through those color-octet channels which can be at the leading 
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order of a s , i.e., the channel with the quantum number 1 So and 3 Pj with J = 0,2. There is a velocity- 
power counting rule to determine the relative importance of NRQCD matrix elements|12( |2] for a given 
quarkonium. In the factorized form the second term is for the relativistic correction of the channel 1 5o- 
The coefficients in the front of the NRQCD matrix elements can be calculated with perturbative QCD, 
they are series in a s , e.g., 

G^CSo) = CfffSo) + GffSo) + • • • > ( 4 ) 

where the subscriber 0(1) stand for the tree(one-loop) contribution. If the factorization holds, these 
perturbative coefficients should not contain any I.R. divergence. To determine these coefficients, one 
replaces the quarkonium H with a QQ state and calculates T^ v and the NRQCD matrix elements. By 
comparing both sides of Eq.(3) calculated with the QQ pair the perturbative coefficients can be extracted. 
If the factorization holds, soft divergences in T^ v will have the same form as those appearing in the matrix 
elements so that the perturbative coefficients do not contain any soft divergence. 

To study the factorization we need to calculate the tensor T^ v in perturbative theory after replacing 
the quarkonium H with those QQ states: 

T^(P,q„QQ) = W d*xe^( |J-(0)|Q(^)Q(^) + X){Q(p 2 )Q( Pl ) +X\J»(x)\0). (5) 
x J 

It should be noted that the heavy quarks carry different momenta in the amplitude and its complex 
conjugated. This will enable us to identify different states of the QQ pair. These momenta are given as: 

V1 = \P + Q, P2 = \P-q, P 'i = \p + q', p' 2 = \p-q'. (6) 
In the rest frame of the QQ, q^ = (0, q) and q 1 ^ = (0, q') with q 2 = q' 2 and 

P 2 = (pi+p 2 ) 2 = 4(m 2 +q 2 ), v=-i v' = ^, (7) 

m m 

v(v') is the velocity of the heavy quark in the rest frame. At tree level, the unobserved state X contains 
only one gluon. By expanding v and v' and identifying quantum numbers, one can determine four 
coefficients in Eq.(3), i.e., 

Tg"(P,q 7 ,QQ) = F^fSo) ■ {Q\O 8 CS , QQ)\0) + G^So) " (^(^o, QQ)\0) 

+F^CP ) • (0|O 8 ( 3 P ,Q0)|0) + F ^( 3 P 2 ) • (0\O^ } ' {kl} eP 2 ,QQ)\0) 
+0(v 4 ). (8) 

In the expansion in v and v', the leading terms give contributions to Fq^^Sq) for the 1 Sq state, the 
next-to-leading terms, which are linear in v and v' like v • v', give contributions to Fq U ( 3 Pq) for the 
3 Po state and Fq V '^ i^P 2 ) for the 3 P 2 state. The next-to-next-to-leading terms are proportional either 
to the tensor v % v 3 or to v n v' 3 . One can decompose the tensor v l v ] or to v h v /:) into the component of 
the S-w&ve with I = and of the D-wave with I = 2. The S*-wave component corresponds to the 
relativistic correction of the 1 5o state. At tree level all these coefficients are nonzero and contain no 
soft divergence, their detailed forms are not important for our purpose because we will show that the 
soft-divergent correction at one- loop to these coefficients is proportional to the tree- level result Tq U . We 
will use Feynman gauge in this letter. In this gauge one can clearly see how soft divergences are cancelled 
or matched in a diagram- by-diagram manner. 
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One-loop corrections consist of two parts. One is the virtual correction, another is the real correction 
in which the unobserved state X contains two gluons or a light quark pair. Beside corrections from 
wave-function renormalization there are many Feynman diagrams. Since we are only interested in soft 
divergences, we do not need to consider all diagrams but only those containing soft divergences. The 
diagrams with soft divergences are given in Fig.l.. To obtain the soft-divergent parts of these diagrams 
we employ the eikonal approximation with some modification. 




Figure 1: Diagrams for contributions containing soft divergences. 

To illustrate the eikonal approximation used here, we consider the contribution from Fig. la to the 
matrix element BP = {Q{p 2 )Q{pi),G{k,e*,a)\J' J '{x)\0}: 
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(9) 



with k 2 = k — k\. The soft divergence appears when ki or k 2 becomes soft, i.e., all components of ki or 
k 2 becomes small, and when ki, hence also k 2 , is collinear to the momentum k of the outgoing gluon. If 
the gluon with ki is soft, the standard approximation is to neglect all ki in nominators and keep only 
the leading order in k\ in the denominators. Therefore the contribution from the soft region of ki can be 
written as: 
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(10) 



We use dimensional regularization with d = 4 — e to regularize divergences. The pole with l/e/ is for I.R. 
divergence with ej = 4 — d. Other poles without the subscriber / are for U.V. divergences. Calculating 
the integral we find that the soft amplitude contains double pole in the form e _1 (e _1 — e^ 1 )' indicating 
that the approximation may not be convenient here. A convenient approximation we will use is to keep 
denominators exact. With this approximation, the soft- and collinear-divergent part of the contribution 
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can be written as 
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and i?j* a — i?^ a sc contains no any soft divergence. In the above equation we have taken a frame in which 

^ = (0, AT, 0,0) = ATn", # = (g+, g -,0,0). (12) 

It should be noted that i? 2(l sc can be written in a covariant form. Performing a similar analysis for other 
diagrams and loop-momentum integration we obtain the soft-divergent part of the one-loop correction to 
The sum of contributions from Fig. la, Fig. Id and Fig.le can be written in a compact form: 



^la,sc + Rld,sc + ^le 
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i.e., the soft-divergent part is proportional to the tree- level amplitude. The double pole in ej is from the 
overlap of collinear- and soft region of the loop momentum. The contribution from Fig. lb and Fig.le 
contain not only soft divergences but also Coulomb singularities. With our approximation the sum of 
these two diagrams gives: 
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where • • • stand for terms which are finite with e, v — > 0, and higher orders in -u. 

Putting everything together, we obtain the soft divergent part of the virtual one-loop contribution to 
T^ v in Feynman gauge: 
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(15) 



where • • • denote finite parts and the correction terms in the third line come from wave-function renor- 
mailization of the heavy quarks and the gluon. 
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Now we turn to the real correction. The real correction consists of amplitudes of two gluons or a 
light quark pair. The two gluons can be emitted by heavy quarks and by gluon splitting, the light quark 
pair can only be generated through gluon splitting. It should be noted that a light quark pair can also 
be produced in other ways, but the contribution of this is irrelevant for the factorization discussed here, 
because in this case the heavy quark pair is in states other than those given in Eq.(8) and the contribution 
at order a 2 s is free from soft divergences. The Feynman diagrams for the real correction is given in Fig. 2. 




Figure 2: Real correction 



We write the amplitude with two gluons in the final state as 
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(16) 



where the first term is only for two gluons emitted by heavy quarks and second is for gluon splitting. If 
two gluons are emitted by heavy quarks, there are only I.R. divergences, when one of the gluons becomes 
soft. Collinear singularities will not appear here because of the massive quarks. When the gluon with 
the momentum k\ becomes soft, the soft contribution of R^ g q with the standard eikonal approximation 
can be written: 
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(17) 



where we only give those contributions explicitly, which are from the color-octet QQ pair with the 
quantum numbers 1 Sq and 3 Pj with J = 0, 2. The contributions from states different than those above 
are indicated by • • •. They are not important for our purpose, but they are important for the complete 
cancellation of soft divergences in the KLN theorem. We will discuss this later. The contribution from 
Rg g q s to T^ v can be calculated easily, we obtain: 
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where we have only collected the contributions from relevant QQ states. Contributions from other states 
and higher orders of v are represented by 
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The amplitude from the gluon splitting is given by: 
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with A; = k\ + /c 2 . The contribution from R% n to T^ v is more difficult to evaluate than that from R^ n n , 

yy^y yy^ 

because there is a overlap between the collinear- and soft region of gluon momenta. However there is 
a standard method, called phase-space slicing method, discussed in detail in |T3[2|. We will use the 
method. If the gluon with k\ is soft, the amplitude can be approximated by: 
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(20) 



This amplitude R gg ^ g ^ s interfered with the R^ g q s will give contributions with the soft divergences which 
are exactly those in the contributions from the interference between R gg q and R gg)9 - By the method 
mentioned above we have the contribution from the interference: 
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(21) 



where s m i n is a cut-off in the phase-space slicing method. Our final result will not depend on it. The 
contributions only from the gluon splitting into two gluons and a light quark pair can be calculated 
directly with the phase-space slicing method. They contain only collinear singularities. The result is: 
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where the same dependence on the cut-off s m i n appears and it will be cancelled by that in T^ u int 
everything together we obtain the infrared divergent part of the real correction as: 
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Finally, we obtain the total soft-divergent part of T^ v at one loop: 
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(23) 



(24) 



where we only give the relevant part in detail. It should be noted that the [■ • •] of the part does not 
contain v • v' from our calculation at the orders considered here. The • • • denote contributions from other 
states of the QQ pair, which can not be produced at tree-level. These contributions are from R gg q, 
represented by • • • in Eq.(18), they contain terms like v • v' because the QQ pair at one-loop order can 
be in a color-singlet P-wave states. 
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Before we turn to our result of NRQCD matrix elements to finalize NRQCD factorization, we briefly 
discuss here the cancellation of soft divergences in KLN theorem. It should be noted that in T^ v defined 
in Eq.(5) the heavy quark pair can be in an arbitrary state which is allowed to be produced at a given 
order of a s , i.e., we do not sum over all possible states of the heavy quark pair. If we sum all these 
states and take v = v', the sum should not contain any soft divergence, as stated by KLN theorem. 
The sum contains not only those terms explicitly given in Eq.(24), but also those represented by • • •, i.e., 
the contributions from those states which can not be produced at tree level, these states are produced 
through emission of two gluons by heavy quarks, whose contributions are represented in Eq.(18) by 
Taking these contributions into account we have the soft-divergent part of the sum: 
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where the second term comes from the contributions of those states which are not given explicitly in 
Eq.(18). By summing of these states and some manipulation the sum can be written into a compact 
form as given in the above. We see clearly that the sum is free from any soft divergence as KLN theorem 
states. 






(b) 



Figure 3: Examples of Feynman diagrams for (a) vertex correction, (b) wave-function renormalization 
and (c) gluon-exchange correction, the broken line is the cut. 

The one-loop correction of those NRQCD matrix elements in Eq.(8) and Eq.(12) can be divided into 
two parts as a virtual part and a real part. The virtual part consists of corrections from vertex and wave- 
function renormailization. Examples of Feynman diagrams at one-loop level for the real part, vertex- 
and wave-function renormailization correction are given in Fig. 3. up to the orders of v we consider. We 
calculate these corrections in Feynman gauge. The one loop correction can be written as: 



(0|O^|0)i = (0|OQ«|0) 



■^■h l+0{v2))+0{v2] 



(26) 



where O®® is any of Og^So), Ps( 1 Sq), 0( 3 Pj) with J = 0,1,2, the matrix element with the subscriber 
0(1) is the tree-level(one-loop) result of the operator. The first term represents the Coulomb singularity, 
the second term starts at the relative order of v 2 , combined with (0|0^|0)o they should be written in the 
form as matrix elements of operators. Using this result we can clearly see that the Coulomb singularity 
in Eq.(24) is reproduced. Hence, one can conclude that all perturbative coefficients F's in Eq.(3) are free 
from the Coulomb singularity and also from I.R. divergences. 

To assess if G^^Sq) contains any I.R. divergence, we need to take the operator mixing between 
Os( 1 Sq) an d others into account which is contained in the second term in Eq.(26). At one loop and 
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order v 2 , the operator can be mixed with the color-singlet operator with quantum numbers 1 P\ and the 
color-octet operator Pg( 1 5'o, QQ). The later is relevant for our case. In Feynman gauge we have 



{0\O 8 (}S ,QQ)\0) 
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(27) 



where contributions at orders higher than v 2 or from mixing of other irrelevant operators and those 
with the Coulomb singularity are represented by • • •. In the above the first line comes from the vertex 
correction, the second from the external quark legs and the third comes from the real gluon exchange. 
With this result, one can see clearly how the soft divergences are absorbed into the matrix element at 
leading order of v on a diagram-by-diagram basis. At order of v° all I.R. divergences are cancelled. At 
the next-to-leading order of v, the net divergence in the above equation exactly matches that in and 
the matching is also on a diagram- by-diagram basis. 





(6) 



Figure 4: Types of Feynman diagrams for contributions with the gauge link. The double line stands for 
the gauge link, (a): Virtual correction, (b): Real correction. 

Our results clearly show that without the modification of relevant NRQCD matrix elements in our 
case the NRQCD factorization holds at order of v 2 and the relativistic correction can also be factorized. 
Adding a gauge link in the those color-octet NRQCD matrix elements, as suggested in the modified 
matrix elements will receive extra contributions. These contributions relevant for our study are given 
by types of Feynman diagrams in Fig. 4. We use Feynman gauge to study these contributions. In other 
gauges, like Coulomb gauge, a gluon can also be exchanged between gauge links. By an infrared power 
counting each contribution from Fig. 4. is with infrared divergences. The virtual contribution from the 
type(a) of diagrams in Fig. 4, after integrating of the energy of the virtual gluon and neglecting terms 
which generate power divergences, is proportional to the integral: 

f d 3 q 1 1 
J (2vr) 3 2|q| |q| + q 3 

In the above the factor l/(|q| + q 3 ) comes from the eikonal propagator of the gauge link, which is 
determined by the moving direction of the quarkonium. We take the direction in the z-direction. It is 



1 + V 5 



|q| + v • q 
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v • q 



(28) 
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clearly that the integral is infrared-divergent. But the real contribution from the type(b) of diagrams in 
Fig. 4 is also proportional to this integral and the proportional coefficient has different sign than that of 
the virtual contribution. The total contribution from the gauge link to the color-octet matrix elements 
is zero at one loop and at any order of v. Therefore, at one-loop level, the suggested gauge link will not 
lead to a violation of the universality of color-octet matrix elements and it will not affect all existing 
one- loop results. 

To summarize: We have examined in detail the NRQCD factorization in inclusive production of a 
quarkonium through e + e _ -annihilation. Our results show that the factorization can be made for pro- 
duction of a QQ pair in color octet with the quantum number 1 So and 3 -Po,2 and also for relativistic 
correction to the 5-wave state. The modification of color-octet NRQCD matrix elements with the sug- 
gested gauge link^J will not affect the NRQCD factorization at one loop in cases studied before and the 
case studied here, because the gauge link does not lead to nonzero contributions to color-octet NRQCD 
matrix elements at one loop. Because of this the universality of these matrix elements holds at one loop. 
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